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Abstract- In this paper, the Laplace Decomposition 
Algorithm (LDA) is introduced for the exact solution of 
the system of weakly singular Volterra integral 
equations. The technique is described and illustrated 
with some numerical applications. The results assert that 
this scheme is rapidly convergent and give the exact 
result using only few terms of its iteration scheme. 
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1, INTRODUCTION 


System of weakly singular Volterra integral 
equations appear in many branches of scientific 
fields[1-6],such as electron emission, radar ranging, 
radio astronomy, atomic scattering, plasma diagnostics 
and optical fiber evaluation, Studies of systems of 
singular integral equations have attracted much concern 
in applied science, 

The system of weakly singular Volterra integral 
equations of the convolution type in two unknowns [7] is 
given by 


U(x) = fx) +f. Tk, QU) + kV (i) ]dt, 


U(x) = £0) +f) [ki QU) + kOx,0V(t)]dt 

(1) 
The kernels k,(x,t), 1 <i, }<2 and the functions f(x), i=l, 
2 are given real-valued functions. The kernel k, are 
singular kernels given by 
k, = 1(x-t)", 1< 1,j< 2. * (2) 
Recall that the kernel is called weakly singular as the 
singularity may be transformed away by a change of 
variable[2]. 


In this work, we used Laplace decomposition 
algorithm(LDA) because this scheme provides the 
solution in a rapidly convergent series with components 
that are elegantly computed. The Laplace 
decomposition algorithm(LDA) was first proposed by 
Khuri[8-9] which is further used by Yusufoglu[10] to 
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solve Duffing equation and Elgazery[11] for Falkner- 
Skan equation. Wazwaz[12] established a necessary 
condition that is essentially needed to ensure the 
appearance of “noise terms” in the inhomogeneous 
equations, The necessary condition for the “noise terms” 
to appear in the components u0 and ul is that the exact 
solution must appear as the part of u0 among other terms. 
The restriction and improvements of Laplace 
decomposition method was given by Khan and 
Gondal[{13]. Zafar et.al, used Laplace decomposition 
method to solve Burger's equaton[ 14]. 


It is worth mentioning that the proposed method is an 
elegant combination of Laplace transform and 
decomposition algorithm. The advantage of this 
proposed method is its capability of combining two 
powerful methods for obtaining exact solution. The aim 
of this work is to establish exact solutions for the system 
of weakly singular Volterra integral equations without 
large computational work. 


Il, THE LAPLACE DECOMPOS 
ALGORITHM 


TION 


In this section, we present Laplace decomposition 
algorithm(LDA) for solving the system of weakly 
singular Volterra integral equations given by (1) and (2). 
The method consists of first applying the Laplace 
transform to both sides of equations in system (1), we 
have 

L{U(x)} =Li{ f(x)} + Lif? [k,,(x.QU(t) + 

k.(x.0V(t) dt}, 

LEV(x)} = Lif (x)} + LYS [ky (xtyU(t) + 
k.Cx.tyVv(t)jdt} (3) 
Using Convolution theorem of the Laplace transform, 
we haveL {U(x)}= L{f(x)}+ L{k,,(x)}L{U(x)}+ 

Lik, GQ) ILI VOX}, 

L{V(x)}=Lif(x)}+ Lik, (x)}L{U(x)} + 

L{k..(x)} L{V(x)} (4) 
Operating inverse Laplace transform on both sides of 
(4), we have 


U(x) = f(x) + L’{L {ky QS) FL {UCR} f+ 

LYLE (x)}L{VOX)} }, 

V(x) = f(x) + EYL {k,,(x) }L{UO)} } + 

LY LAK OL {VON} (5) 
The Laplace decomposition algorithm(LDA) defines 
the solutions U(x) and V(x) by the infinite series 

U(x) = E5..U,, and V(x) = Div, (6) 
Substituting (6) in (5), we have (7) 

Deo U, = f(x) + LAL fk, (X)}L{ E50 U, 3} 


HLL KO) LA Zio Vas 

Tho V, = 80) + LL {KO} L{ 

+ LYUL{KAX)JL{Zeo Va} () 
In general, the recursive relations are given by: 
U,(x) = f(x) and V,(x) = H0x) (8) 


Unis) = LAL EK (x)PL{ EE,U, f+ 
L'{L{k,(x)}L{ D",V,}}, n20 and 
Voor (x) = L'{L (ka (XL A Zi0U, H+ 
L'{L{k,.(x)}L{ ,V,}} n=O. (9) 


Ill, APPLICATIONS 
In this section, some applications are given in order to 
demonstrate the effectiveness of Laplace decomposition 
algorithm (LDA) to solve the system of weakly singular 
Volterra integral equations. 
A, APPLICATION: | 
Consider the system of weakly singular Volterra integral 
equations 
U(x) = bx? = 5x72 +f TU (t(x-1) 
+ V(t)/(x-t)! "Jt, 
V(x) = 1 = x0 = Sx°/2 + fPU(U/(x-t)!* + 
V(ty/(x-t)!"Jdt (10) 
Applying the Laplace transform to both sides of (10), we 


have < 
L{U(x)}= I/s + 2/s’ — 20(4/S)/s° + 
LEf [UQV/Cx-t)" + V(D/(x-t)" Tat}, 
L{V(x)}= M/s - 2/s’= 20(4/5)/s"*+ 
Lif TUC/x-0)'" + V(/x-t)"Jdty (yy 
Using Convolution theorem of the Laplace transform, 
we have 

L{U(x)} = 1/s + 2/s’ — 20(4/5)/s”* + . 

L{U(x)} + P(4/5)/s* L{V(x)}, 

1/s - 2/s'— 2T(4/5)/s”* + 

L{U(x)} + P(4/5y/sL{V(x)} (12) 
Operating inverse Laplace transform on both sides of 
(12), we have 

U(x)=1 +x" -5x"/2 +L {1 (4/5)/s"°L {U0} } + 
LUP(4/5y/s°L{V(x)} f, 

V(x) = 1 =x" = 5x"/2 +L" {P(4/S)/s°L {U(x)} } + 
L'(T(4/5)/sLAV(x)} 


(13) 


The Laplace decomposition algorithm (LDA) defines 
the solutions U(x) and V(x) by the infinite series 


U(x) =E5, U, and V(x) =Ziv, (14) 
Substituting (14) in (13), we have 

Eh) U,=1t x - 5x2 + 

L’ {T(4/Sy/s"L 45%, U, }}+ L" 

{P(4/5\/s°L4 Teo V, }}s (15) 


By V=1-x'- 5x72 4+ L 
{T(4/5V/s"L{ Ero, }}+ L" 
{T(4/5\/sL{ TEV, }} 


From (15) our required recursive relations are given by; 


U,(x)=1+x°-5x""/2 and V,(x)=1-x"-5x"'/2 (16) 
U,, 00 = LED GSy/s"L{U,}} +L" 
(1(4/S)/s"L{V,}}, 120 and (17) 


Vu (x) = L{r(4/sy/s“L{U,}} +L" 
{P(4/S)/s“L{V,}}, n20. 

The first few components of U,(x) and V,(x) by using 
recursive relation (17) as follows immediately. 

U(x) = 5x""/2 — 25[P(4/5)]'x° /6(3/5) and 


V(x) = 5x"°/2 ~ 25[F(4/5) x" /61(3/5) (18) 
Now, by noise terms phenomena, we have 
U(x) = 1 +x? and V(x) = 1 - x? (19) 


Which are the exact solution and satisfy the (10). 

B. APPLICATION: 2 

Consider the system of weakly singular Volterra integral 
equations 

U(x) = x +x) = 25x77 +f [U(/(x-t)"* 4 

V(t)/(x-t)'"Jdt, 

V(x) =x =x = 25x12 +f [U/(xt) + 

V(tV(x-t) Jat (20) 


Applying the Laplace transform to both sides of (20), we 
have 


L{U(x)} = I/s? + 2/s? - 210°(2/5)/s"* + 


vs 


Life (Uix-ty* + V(t)/(x-t)* Jade}. 
L{V(x)} = l/s’ - 2/s’ - 27°(3/5)/s'°+ 
Lif" [U(Dx-t)* + V(t)/(x-t) dt} (21) 
Using Convolution theorem of the Laplace transform, 
we have 

L{U(x)} = L/s* + 2/s° - 20(2/5)/ 
T(2/5V/s*L{U(x)} + 1(2/SV/s“L{V(x)}, ( 
L{V(x)} = I/s" - 2/s’ - 21F(3/S)/s"* + 
T(3/5)/s““L{U(x)} + T3/S\/s“L{V(x)} 

Operating inverse Laplace transform on both sides of 
(22), we have 

U(x)=x +x -25x7/7+ LT (2/5)/s* 
L{U(x)}} +L {P(2/Sy/s* L{V(x)}}, 
V(x)=x- x -25x"/12 +L" {1(3/5)/s° 
L{U(x)} } + LTG/SYs"L{VOX)}} 


re) 
Nn 


(23) 


25 


The Laplace decomposition algorithm(LDA) defines 
the solutions U(x) and V(x) by the infinite series 


U(x) =ZTiq U, and V(x) =E% ov, (24) 
Substituting (24) in (23), we have 

Thy U=x tx? - 25x"/7 + L" {T(2/5)/s™* 

L{ZzoU, }}HF LL (2/S)/s"L { Zz V, }}, (25) 


Dry VHX x = 25x12 + LY (T(3/S)/sL 
(27,U, }}4 L" {P(3/5)\/s°LEEA LV, $4 


From (25) ourrequired recursive relations are given by: 
U,Q)=x+x°-25x"°/7 and 


V(x) =x— x? - 25x°/12 (26) 
U,.00 = L{P(2/5)/s*L{U,}} + 
L'{(4/5y/sL4V,}}, n20 and 

V,.(x) = L'{D(4/Sy/s“L fU,} } + 

L'{r(4/sy/s’L fV,} }, n20. (27) 


The first few components of U,(x) and V,(x) by using 
recursive relation (27) as follows immediately 

U(x) = 25x"°/7 - 25(P(2/5)}°x"*/181 (4/5) 
-1(2/5)P-(3/5)x° and 

V(x) = 25x12 - 125[7(3/5)]x"/33F 


(1/5) - PQ/SEB/S)x° (28) 
Now, by noise terms phenomena, we have 
U(x)=x+x° and V(x)=x-x* (29) 


which are the exact solution and satisfy the (20). 


IV. CONCLUSION 

In this paper, we have successfully developed the 
Laplace decomposition algorithm (LDA) for the 
solution of the system of weakly singular Volterra 
integral equations, The given applications showed that 
the exact solution have been obtained even with just the 
two first terms of the LDA solution, which indicates that 
the proposed method LDA need much less 
computational work. The proposed scheme can be 
applied for the system more than two weakly singular 
Volterra integral equations, 
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